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Abstract. It has been suggested that relativistic jets in quasars may contain a considerable amount of thermal matter. In this 
paper, we explore the possibility that the Ka line from the thermal matter may appear at tens of keV due to a high Doppler 
blue-shift. In the jet comoving frame, the energy density of photons originally emitted by the accretion disk and reflected off 
the broad line region clouds dominates over that of photons of other origin. We discuss the photoionization states of the thermal 
matter and find that the irons elements are neutral. The high metallicity in quasars enhances the possibility to detect the thermal 
matter in the relativistic jet in some radio-loud quasars. A highly Doppler blue-shifted Ka line may be detected. We make a 
prediction for 3C 273, in which the Ka line luminosity might be of the order 3.0 x 10 44 erg s with an equivalent width of 
2.4 keV. Such a line could be detected in a future mission. 

Key words, galaxy: radio galaxies: active - quasars - individual: 3C 273 



C^* 1. Introduction 



It is well known that the observed powerful relativistic jets are 
an important ingredient in radio-loud quasars. The non-thermal 
electrons are responsible for the multiwaveband radiation al- 
though the specific acceleration mechanism remains a matter 
of debate. Little attention has been paid to questions related to 
the thermal matter in the relativistic jets. 

As an example of a galactic jet source, SS433 has been ex- 
tensively studied both observationally and theoretically (Kotani 
et al. 1997; Brinkmann & Kawai 2000). Chandra discovered 
several lines from highly ionized atoms, such as Fe XXV, Fe 
xxiv, Co xiv, S xvi, Lya and Ly/3, Ne x and Mg xi etc. 
(Marshall et al. 2002). The chemical composition in these jets 
definitely includes heavy elements rather than pure pair plasma. 
This could be explained by a model in which the emission 
lines are originally from the hot plasma expanding in the jet 
(Brinkmann & Kawai 2000; Memola et al. 2002). In radio-loud 
quasars, the case is highly uncertain. A heavy jet mainly com- 
posed of proton-electron plasma has been suggested by Celotti 
& Fabian (1993) based on the kinetic luminosities of jets found 
in a large VLBI sample. The absences of soft X-ray bumps in 
radio-loud quasars lead to the exclusions of a pure pair and pure 
proton-electron plasma, most likely, the relativistic jet is pair- 
dominated numberwise but still dynamically dominated by pro- 
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tons (Sikora & Madejski 2000). The measurement of polariza- 
tion in a few objects seem to favor the pair plasma (Wardle 
et al. 1999, Hirotani et al. 1999), but the linear polarization 
strongly supports a normal plasma as the main composition 
(Fraix-Burnet 2002). Ruszkowski & Begelman (2002) find that 
the electron-proton and electron-positron jets can lead to the 
same circular and linear polarization in 3C 279. 

There is growing interest in the presence of thermal matter 
in relativistic jets (Celotti et al. 1998). The thermal matter, as 
argued by Celotti et al. (1998), may be due to: 1) not 100% 
matter can be accelerated to relativistic energy; 2) non-thermal 
matter cools down and is thermalized before being reacceler- 
ated; 3) some thermal matter might also be trapped at the base 
of the jet as they form and some are loaded by the surround- 
ing external medium (but see Lyutikov & Blandford 2002 for 
a different view). Kuncic et al. (1997) used CLOUDY to make 
detailed calculations of emission lines from the thermal matter 
immersed in the non-thermal radiation field in the jet. The ba- 
sic features are the presence of emission lines in the extreme 
ultraviolet band. However, the situation in radio-loud quasars 
may not be so simple. There are three types of possible sources 
for the ionizing photons in the relativistic jet: 1) synchrotron 
photons (Blandford & Konigl 1979); 2) accretion disk photons 
(Melia & Konigl 1989; Dermer & Schlickeiser 1993); 3) dif- 
fuse photons in the broad line region (Sikora et al. 1994). At 
present the observations do not allow to decide which might 
dominate in blazars. In the jets comoving frame the thermal 
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matter sees photons originating from the disk and reflected by 
the clouds in the broad line region. These photons may dom- 
inate over the local synchrotron photons (Sikora et al. 1994). 
Moreover, it is of great interest to note that a high metallicity is 
common in quasars (Hamann & Ferland 1999). 

In this paper, we show that the thermal matter in the jet 
is neutral and the observational features of the thermal matter 
will mainly be the presence of a highly Doppler shifted Ka line, 
which could be detected in some radio-loud quasars by future 
instruments. 

2. Emission from Thermal Matter in Jets 

2. 1. Ionization of thermal matter and Fe Ka line 

Detailed constraints on temperature, density and size of the 
thermal matter in relativistic jets have been discussed by Celotti 
et al. (1998), who suggest that the thermal matter may exist 
as cold clouds. The microscopic properties of the surviving 
clouds in the relativistic jet depend on many processes. The 
density and size of the clouds are largely uncertain. In this pa- 
per, we take typical values for the density, temperature and 
the size as given by Celotti et al. (1998), e.g. n, = 10 14 n U 4 
cm 3 , T = 10 5 r 5 K and R th = l0 9 R th , 9 cm. For illustration, 
we consider the clouds in the jet at a distance R from the 
central black hole, with R = 10R g = 1.5 x 10 15 cm, where 
R g = 1.5 x 10 14 (cm)m9 and mg = M B h/1O 9 M . We approx- 
imate the Doppler factor D with the Lorentz factor F and take 

r 10 = r/io*i. 

The peak frequency of the thermal emission from such 
clouds in the jet will be TkT w 0.1 (keV)Fio75 in the ob- 
server's frame, where k is Boltzmann constant. If such a soft 
X-ray emission is absent in the observed continuum, the num- 
ber of the cold clouds No should be constrained by 

N < 5.62 x 10 8 L S x,46r m 4 r 5 4 ^ 9 . (1) 

where the soft X-ray luminosity Lsx = 10 46 Lsx,46 ergs -1 . 
Assuming the jet geometry to be a cone with opening angle 
6j ~ 1/r, the cross section radius of the jet is a = 0jR . 
Since the sound velocity inside the cold clouds is much less 
than that in the intercloud medium, the expansion of the clouds 
can be neglected. The filling factor of the thermal clouds in 
the jet volume is / « £f° (4n/3)R 3 fh / {nR 3 /3r 2 ) * 6.66 x 
1(T 8 where R l5 = fl/10 15 cm and the opening an- 

gle is taken as 6 « 1 /T, indicating that the jet is very clumped 
with thermal clouds as suggested by Celotti et al. (1998). These 
general constraints on the thermal matter are in agreement with 
those of Sikora et al. (1997). 

According to the disk-corona model by Haardt & Maraschi 
(1993), most of the gravitational energy will be released in the 
hot corona as X-ray emission, and UV emission as reprocessed 
X-rays. We take a similar scenario for the accretion disks in 
flat spectrum radio quasars. As an example, we take the accre- 
tion rate m = Lboi/^Edd = 0.2, where the Eddington luminosity 
is ^Edd = 1-26 x 10 47 m9(erg s _1 ). The bolometric luminosity 
is then Lboi = 2.52 x 10 46 (erg s _1 ). For simplicity we assume 
that half of the bolometric luminosity is released in X-rays, 



L x = 1.26 x 10 46 m 9 (erg s" 1 ), and the other half in optical/UV, 
Luv = 1.26 x 10 46 m9(ergs s _1 ). Some of the disk emission will 
be reflected by the clouds in the broad line region. 

Sikora et al. (1994) point out that the photons in the rel- 
ativistic jet received from diffuse scattering of accretion disk 
photons in the BLR may dominate over the local synchrotron 
photons and the photons received directly from the disk. The 
nature of the BLR clouds remains open (Alexander & Netzer 
1994; Baldwin et al. 2003), the covering factor £ w 0.1 is 
indicated by the energy budget of emission lines from BLR 
clouds (Netzer 1990). The typical density and temperature are 
n e = 10 n ~ 12 cm 3 and T = 10 4 ~ 5 K in the clouds of BLR, 
and the ionized fraction of dimension is typically ( = 10 12 cm. 
The Thomson scattering optical depth is r es = n e cr T { ~ 0.1, 
where crj is the Thomson cross section. Collin-Souffrin et al. 
(1996) calculate emergent spectra from an optically thin cloud 
(with solar abundances) radiatively heated in detail. They found 
that about 1 % of the incident radiation will be reflected for a 
cloud with column density "L = 10 22 g cirT 2 (Thomson depth 
r es w 0.006) and the reflected amount is insensitive to the 
cloud's temperature. In the optically thin regime, the reflected 
flux will be approximately proportional to the optical depth. 
Therefore, the reflected fraction of the incident radiation from 
disk is "R > 0.1 for the BLR clouds with r es -0.1. The en- 
ergy density of the diffuse X-rays is Ux = L 1 ^ /47i7? 2 LR c, where 
L^ f = ^KL X , Rblr is the scale of the BLR and c the light speed. 
In the jet frame, the energy density of the received X-rays by 
the thermal matter is then 

U' x = F 2 U X = 0.9 (erg cm^F 2 ^/^, (2) 

where we take the reflected fraction of hard X-rays l*R = 0.01 
conservatively, the BLR size ^o.ipc = ^blr/0.1pc and L^ f 44 = 
L^ f /10 44 erg s~'. The ionization parameter defined by 

3 = ^ - 6.0 x 10- 4 T\ a Ll%R^ w n;\ A T-\ (3) 

indicating that most of the atoms are neutral (Krolik & Kallman 
1984) and the iron Ka line is at 6.4keV. It is thus expected that 
the thermal clouds will maintain their temperatures. 

2.2. Iron abundance and line luminosity 

As shown by Monte Carlo simulation of Reynolds (1996), the 
iron Ka line is the most prominent among the emergent lines 
in the reprocessed spectrum. We neglect the lines of other el- 
ements since they are much fainter than the iron Ka line. We 
take the relative abundance of iron as J[ = Zqso/Z , where 
Zqso is the iron abundance of QSO and Z = 3 x 10~ 5 . Krolik 
& Kallman (1987) (see also Liedahl 1999 for details) show that 
the Fe K-edge opacity is 

kko = 0.67* es /(E)W£-J (4) 

where K es = 0.34 is the Thomson scattering opacity, and /(H) 
is a slow function of H. The optical depth of one cold cloud is 

t k = 0.45/(E)^ 10 « t ,i4flt,9, (5) 
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where JI\q = Ji/10. For simplicity we assume that the reflected 
spectrum as the photoionization source is in a simple form as 
L' E = L' (E/E Q y y within E Q < E < E u where y « 0.7 is 
the index of the reflected spectrum by the clouds in the BLR. 
In the jet's comoving frame, J L' E dE = T 2 L^ f , yielding L' () = 

(1 - y)CoE l T 2 Lf(EJE r\ and C = [(Ei/£o) 1_r - 
The reflected spectrum is assumed to be the same as the inci- 
dent radiation field. This may be not exact, but the accuracy 
is enough for the present goal since the Ka line flux mainly 
depends on the total energy of hard X-rays. 

With the optical depth given by equation (5), the luminosity 
of the Ka line from A^o clouds in the jets comoving frame will 
be given by 



L Ka *(Y)T K L Q —(—j 



(Nq6lo\ 



(6) 



where (Y) is the fluorescent yield for the production of Ka, 5a> 
is the solid angle of one cloud subtended at the X-ray contin- 
uum source. In the observer's frame, the received luminosity is 

Igj = 3.0 x 10 44 (erg/s)rf J?liont,i4<, 9 ^r5 2 ^44 ( 7 ) 

with (Y) * 0.6, /(H) = 1.7, E Q * E Ka , Ei/E = 100/6.4 and 
A^o = 5.62 x 10 8 . The observed energy of the Ka line will be 
* 64 (keV) T 10 . 
The profile of such a line may be mainly broadened by the 
relative motion among the cold clouds. As argued by Celotti 
et al. (1998), the strength of a comoving magnetic field is of 
typical value B — 2 x 10 3 L^ 46 R^T^ Gauss, where Lj e t,46 = 
^jet/10 46 erg s _1 is the power of the jet as Poynting flux. Such a 
magnetic field can confine the relative motion among the cold 
clouds with respect to the relativistic bulk flow, otherwise the 
collimation of the jet will be broken down. This random veloc- 
ity of the clouds v c can be estimated by ^n,m p v\ = B 2 /8n. 
For the typical value, we have v c /c = 1.5 x 10~ 3 , namely, 
AE « rE Ka v c /c w O.lFiokeV. If the thermal clouds follow 
the opening angle of the jet due to random motion, on the other 
hand, their relative velocity would be of order of v c /c « O.ir^J. 
This causes a broadening of AE « 6.4keV. The resolution of 
the line profile may probe more detail dynamics of the jet in 
future. 

From eq. (7), we see that the observed line luminosity is 
very sensitive to the Lorentz factor F and proportional to the 
iron abundance J[. We use the maximum number of cold clouds 
(see equation 1), the predicted luminosity of the iron Ka line is 
the upper limit. This limit is due to the absence of features of 
bulk Comptonization in the soft X-ray band. We note that the 
above model only works within BLR (~ O.lpc). 

2.3. Candidate objects with a blue-shifted Ka line 

The proposed model naturally relates to the y-ray emission 



model advocated by Sikora et al. (1994). L y /L syn w U^/U B 
and t/yy > t/ syn implies that L y /L syn > 1 if the equipartition 
between magnetic field and relativistic electrons Ub = t/ e is 
fulfilled. Here L y is the y-ray luminosity, L syn is the luminos- 
ity due to synchrotron emission, f/ B and are the energy 



densities of the magnetic field and of the reflected UV lumi- 
nosity, respectively. CGRO observations of y-rays from blazars 
show the importance of reflection (Dondi & Ghisellini 1995). 
The potential candidates for a blue-shifted iron Ka line should 
be those objects in which L y /L syn > 1. On the other hand the 
present model needs metal-rich thermal matter. This condition 
corresponds to the observable indicator of metallicity as N v/C 
IV> 0.4 (Hamann & Ferland 1999). We thus have the criteria 
for potential candidates as 



-^>i; 



N V 
CIV 



> 0.4. 



(8) 



These conditions are likely satisfied in flat spectrum radio 
quasars. The superluminal motion of y-ray bright blazars show 
that T ~ 10 is common among the y-ray loud AGNs (Jorstad 
et al. 2001). In the sample of Dondi & Ghisellini (1995), the 
mean ratio of L y /L op « 30. We expect that the highly Doppler- 
shifted Ka line appears in most of these objects with high flux 
ratio N v/C IV in the sample of Dondi & Ghisellini (1995). 

3. A predicted blue-shifted Ka line in 3C 273 

3C 273 is a typical blazar with strong emission lines, power- 
ful continuum emission from radio to y-rays (see an exten- 
sive review of Courvoisier 1998). A very strong MeV bump 
has been discovered in 3C 273 by the Compton Gamma-ray 
Observatory (Lichti et al. 1995). This MeV feature could be 
explained by several different models, for example, a break in 
the electron injection function (Ghisellini et al. 1996), incom- 
plete cooling of relativistic electrons (Sikora et al. 1994); and 
pair cascade process (Blandford & Levinson 1995). However, 
in such models, one encounters other difficulties to reconcile 
the broad band multi-wavelength continuum emission [see a 
brief review of Sikora et al. (1997)]. As argued by Sikora et 
al. (1997), the MeV feature could be explained naturally by 
the "hot electrons" version of the external radiation Compton 
model (Sikora, Begelman & Rees 1994) provided the plasma 
is not dominated by pairs. The future detection of the highly 
blue-shifted Ka line could help clarify the composition of the 
relativistic jet in radio-loud quasars. 

The mass of the black hole in 3C 273 can be esti- 
mated from the absolute magnitude Mr of its host galaxy via 
log Mbh/Mq = -0.5M R - 2.96, M BH = 1.2 X 10 9 M G from 
M R = -24.4 (McLure & Dunlop 2001). The accretion rate of 
the black hole can be roughly estimated from the big blue bump 
(Walter et al. 1994, Wang & Zhou 1996). We find m * 0.4, 
about half the Eddington luminosity, similar to Courvoisier 
(1998). The metallicity can be estimated through the flux ra- 
tio of N V/C IV (Hamann & Ferland 1999). The observed flux 
ratio of N v/C IV is 0.46 in 3C 273 (Baldwin et al. 1989), 
which gives a metal abundance = Z/Z Q « 10 from Fig. 6 
of Hamann & Ferland (1999). The superluminal motion has 
been extensively studied, the latest observation shows that the 
apparent velocities for different components are from 9c to 22c 
(Jorstad et al. 2001). Here, we take F « 10. 

For the parameters of 3C 273, we estimate a Ka line 
luminosity L* s * 3.0 x 10 44 erg s" 1 and the flux F Ka = 
3.9 x 10~ 12 erg crrT 2 s _1 . For the equivalent width, we use 
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Fig. 1. A plot of the line profiles with future missions. A Gaussian 
profile for the line is assumed in the jet comoving frame. We assume 
two cases with energy resolution of AE = IkeV and AE = 6.4keV, 
respectively. The underlying continuum is the INTEGRAL spectrum 
observed by Courvoisier et al. (2003). 

the INTEGRAL continuum spectrum, F E = 3.6 x 10 "£7 073 
erg ctrT 2 s keV" 1 (Courvoisier et al. 2003). The equivalent 
width is given by EW(Ka) * F Ka /F E * 2.4 (keV). 

Using the observed background level for SPI, both analyt- 
ical calculations as well as Monte Carlo simulations show that 
SPI will not be able to detect such a line in a total observing 
time of one million seconds. Figure 1 shows that for a next 
generation instrument for which we assume a 1 keV energy 
resolution a signal to noise ratio of ~ 1-3 per resolution el- 
ement is necessary in order to detect the line. This represents 
an improvement by a factor of 3 ~ 4 on the presently available 
data that we expect will be available in missions like NeXT 
(New X-ray Telescope, Takahashi et al. 2004) and XEUS (X- 
ray Evolving Universe Spectrometer). The line profile for a 
well-collimated jet with the maximum width (AE = 6.4keV) 
discussed in §2.2 is also plotted in Figure 1. 

4. Conclusions 

We show that future instruments may allow us to detect a ther- 
mal line emitted by matter in a jet directed towards us and 
therefore shifted to the blue by a factor that reflects the gamma 
factor of the jet. This would allow a direct measurement of the 
jet gamma factor and give very important indications on the as 
of yet not clear nature of the jet. The line flux is determined 
by 3 factors: the ionization level, the Lorentz factor and the 
metal abundance. Using realistic values for these parameters, 
we show that the next generation hard X-ray instruments may 
well measure this component. 
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